Quantum spin liquid near Mott transition with fermionized 7r-vortices 
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In this paper, we study the non-magnetic insulator state near Mott transition of 2D 7r-flux Hub- 
bard model on square lattice and find that such non-magnetic insulator state is quantum spin liquid 
state with nodal fermionic excitations - nodal spin liquid (NSL). When there exists small easy-plane 
anisotropic energy, the ground state becomes Zi topological spin liquid (TSL) with full gapped 
excitations. The U(l) x (7(1) mutual-Chern-Simons (MCS) theory is obtained to describe the low 
energy physics of NSL and TSL. 
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I. INTRODUCTION 

People have been looking for quantum spin liquid 
states in frustrated spin models for more than 20 yearsi^. 
For example, various approaches show that quantum spin 
liquids may exist in two-dimensional (2D) S — 1/2 J1-J2 
model or the Hcisenbcrg model on Kagomc lattice^"—. In 
these models, the quantum spin liquids are accessed (in 
principle) by appropriate frustrating interactions. How- 
ever, the nature of the quantum disordered ground state 
is still much debated^—. People have guessed that the 
quantum disordered state may be either algebra spin liq- 
uid stat o 16 ' 17 , or Z2 spin liquid stated or chiral spin 
liquid« 

On the other hand, the experiments in the organic ma- 
terial k- (BEDT - TTF) 2 CU 2 (CN) 3 indicate the real- 
ization of the quantum spin liquid states 2 ^—. Then, 
U(l) and SU(2) slave-roton theories of the Hubbard 
model were formulated on the triangular or honeycomb 
lattice a 24 ' 25 . Recently, the quantum spin liquid state near 
Mott transition of the Hubbard model on honeycomb lat- 
tice has been conformed by different approaches 2 ^—. In 
particular, in Refj 2 ^, quantum spin liquid state has been 
predicted by quantum Monte Carlo (QMC) simulation. 

In Refi2£, the possible non-magnetic state in the nodal 
antifcrromagnetic (AF) insulator (an insulator with AF 
spin density wave ordering and massive Dirac fermionic 
excitations) is predicted to near the Mott transition of 
the 7r-flux Hubbard model (See Fig.l). Such type of 
quantum disordered states in bipartite lattices is also 
not driven by frustrations, as people have done in varied 
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spin models. Instead, they come from quantum fluctu- 
ations of a relatively small effective spin-moments near 
Mott transition. What's the nature of the possible non- 
magnetic state here? In this paper we will answer the 
question. We find that for the isotropic case, the non- 
magnetic state with SU(2) spin rotation symmetry is a 
new type of spin liquid - nodal spin liquid (NSL) with gap- 
less fermionic excitations and roton-like excitations; on 
the other hand, for the anisotropic case by adding small 
easy-plane anisotropic term, the non-magnetic state with 
weakly breaking SU(2) spin rotation symmetry becomes 
a Z2 topological spin liquid (TSL) with topological de- 
generate ground states. 

The paper is organized as follows. In Sec. II, we study 
the spin fluctuations of the nodal AF insulator (NAI) in 
7r-flux Hubbard model on square lattice based on a formu- 
lation by keeping spin rotation symmetry. In Sec. Ill, we 
study the properties of half skyrmions - topological soli- 
tons of the nodal AF insulator and show their induced 
quantum numbers and statistics. In Sev.IV, nodal spin 
liquid is proposed to be the ground state of the quantum 
non-magnetic insulator state in the nodal AF insulator. 
In this section, we use U(l) x U (1) mutual-Chern-Simons 
theory to learn the properties of NSL and TSL. Finally, 
the conclusions arc given in Sec.V. 



II. NON-MAGNETIC STATE IN NODAL AF 
INSULATOR OF THE HUBBARD MODEL ON 
7T-FLUX LATTICE 

In this paper we will focus on the 7r-flux Hubbard 
model on square lattice^. The Hamiltonian of it is 



H = — 



tij C i C 3 



h. 



(1) 



FIG. 1: The scheme of the spin liquid state near Mott tran- 
sition of the 7r-flux Hubbard model 



Here Ci = {ci^,Ci\) T are defined as electronic annihi- 
lation operators. U is the on-site Coulomb repulsion. 

denotes two sites on a nearest-neighbor link, 
and fiii are the number operators of electrons at site i 
with up-spin and down-spin, respectively. There is a n- 
flux phase when a particle hops around a plaquette in 
a 7r-flux lattice. So the nearest-neighbor hopping tij in 
a 7r-flux lattice could be chosen as £4,1+3 = t, t^-s+j) = 
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te ±l %— . Although 7r-flux phase docs not break trans- 
lational symmetry, we may still divide the square lat- 
tice into two sublatticcs, A and B. After transform- 
ing the hopping term into momentum space, we obtain 
£k = ±^/4t 2 (cos 2 k x + cos 2 k y ). So there exist two nodal 
fermi-points at ki = (§,§), ks = (§,—§) and the spec- 
trum of fermions becomes linear in the vicinity of the 
two nodal points. In the non-interacting limit, the 7r-flux 
Hubbard model is reduced into a free fermion model with 
nodal fermions at half filling. 

Because the Hubbard model on bipartite lattices is 
unstable against antifcrromagnetic instability, at half- 
filling, the ground state may be an insulator with AF 
SDW order with increasing interacting strength. Such 
AF SDW order is described by the following mean field 
order parameter ((— l) l c\ir z Ci) = 4^. Here M is the stag- 
gered magnetization. Then in the mean field theory, the 
Hamiltonian of the 2D 7r-flux Hubbard model is obtained 
as 

« = - £(*M4cj + h.c.) - (-I)' *4<*A (2) 

where A = ^M- leads to the energy gap of electrons and 
o z is the Pauli matrix. After diagonalization, the spec- 
trum of the electrons is obtained as = ± + A 2 . 
By minimizing the free energy at zero temperature, the 
self-consistent equation of (|2|) is reduced into j? ^2 JE~ = 

k 

1 where N is the number of the particles. 

From mean field approach, one may see that the MI 
transition of the 7r-flux Hubbard model occurs at a crit- 
ical value about J7/t ~ 3.11. In the weakly coupling 

I 



and the transverse spin susceptibility is 

x± = [( ^?«^bjf r, - 2ur ' (6) 

See detailed calculations in appendix. In the paper we 
will set the spin velocity to be unit. 

The coupling constant g of the 7r-flux Hubbard model 
has been obtained in Ref.— . In Refi^ 3 ., a quantum crit- 
ical point at g c = ^ was found that corresponds to 
U/t ~ 4.26 of the 7r-flux Hubbard model. For the case of 
g > g Cl a quantum non- magnetic insulator (NMI) state 
with gapped spin excitations as m s = 47r(-^ — i) appears 
that corresponds to the region of 3.11 < U/t < 4.26; 
For the case of g < g c , the ground state is a long range 
AF order that corresponds to the region of U/t > 4.26. 
Here the momentum cutoff is introduced as A = min(a , 



limit (U/t < 3.11), the ground state is a semi-metal (SM) 
with nodal fermi-points. In the strong coupling region 
(U/t > 3.11), due to M / 0, the ground state becomes 
an insulator with relativistic massive fermionic excita- 
tions. However, the non-zero value of M only means the 
existence of effective spin moments. It does not necessar- 
ily imply that the ground state of NAI is a long range AF 
order in the mean field theory. Thus one needs to exam- 
ine stability of magnetic order against quantum fluctu- 
ations of effective spin moments based on a formulation 
by keeping spin rotation symmetry, a z — > • a. 

Within the haldane's mapping, the spins are 
parametrized as f2i = (-l) l n iv /l - L 2 + Lf22r— . Here 
n, is the AF order parameter and |nj| = 1, Lj is the 
transverse canting field, which is chosen to Lj • n, = 0. 
By replacing the electronic operators c] and Cj by Grass- 
mann variables c* and Cj , we get the effective Lagrangian 
with spin rotation symmetry as 

£ e ff = £ c*d T c i —y^(ti ! jC*Cj+h.a)-y^ j (~1) % Ac*fti-aci. 

» (ij) » 

(3) 

After integrating the massive fermions and transverse 
canting field, an effective 0(3) nonlinear a model appears 
that describes the long- wavelength spin fluctuations with 
spin rotation symmetry-£ ' 33 i 38 ~ — 

C s = 1- (d^nf . (4) 



The coupling constant is obtained as g = ./ — where 



(5) 

I 

2A) 38 i 39 (a is lattice constant). The existence of a non- 
magnetic insulator state provides an alternative candi- 
date for finding spin liquid state. An interesting issue is 
the nature of the non-magnetic insulator. Is it a valence- 
band crystal, or algebra spin liquid state, or a new type 
of quantum state? In the following parts we will answer 
the question. 

III. HALF-SKYRMION AS FERMIONIC 
EXCITATION 

To learn the nature of the non-magnetic insulator, one 
may study its excitations. In the non-magnetic insulator, 
we get the effective model of massive spin-1 excitations 

£ s = l[(^n) 2 + m 2 n 2 ]- (7) 



Ps = 



N 4-^ 



t 2 [cos (2k x ) (A 2 + 8t 2 + At 2 cos (2k y )) + A 2 + 3t 2 + i 2 cos (Ak x )} 
2(|&| 2 + A 2 )f 
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Or using the CP(1) representation, we have 



io/i)z|" 



2 21 



(8) 



Here is intro- 



where z is a bosonic spinon, z = (zi,Z2), n; 
zz = 1, a M = -|(z5 M z - <9 M zz). 
duced as an assistant gauge field. Because the bosonic 
spinons have mass gap, we may integrate them and get 
£ e ff = ■^■(dfj.a^) 2 with e 2 — m z = |m s . Due to the 
instanton effect, the gauge field a M obtains a mass gap 
and bosonic spinons that couple the gauge field a p are 
confined. Now the lowest energy excitations are gapped 
spin wave. However, the answer is not quite right. Ap- 
plying the Oshikawa's commensurability condition and 
Hastings' theorem to the present case, the quantum state 
without spontaneously symmetry breaking can be either 
a topological order with degenerate ground state or a 
uniform ground state with triplet excitation gap must be 
accompanied with other gapless excitation o 47 i 48 . In this 
paper we indeed find that non-magnetic insulator is ei- 
ther topological spin liquid with degenerate ground state 
or a nodal spin liquid with gapless fermionic excitations. 



the other a down-spin polarized core Q = ^, 

X(x-xq) A(y-yo) 



(n hs ) 2 



(12) 



V|r-r | 2 + A 2 Vlr-roP + A 2 

* ). 

v/|r-r | 2 + A 2 

Fig. 2 and Fig. 3 show the two types of merons. From 
them, one can see that the topological charge of a half- 
skyrmion is determined by both the spin configuration 
and the polarized direction of AF order in the core. 
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FIG. 2: The scheme of (nh s ) 15 the meron with topological 
charge —1/2. up-spins locate at the center. 



In the following parts, we focus on the half-skyrmions 
(topological vortices with half topological charge), 



Q = 



d 2 r^-e Ql/X n-d v nxd x n=±l. (9) 

47T Z 
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In the vector n representation, the solutions of the half- 
skyrmion of the continuum limit are 4 *^— 



n hs = ( 



±- 



X[x - Xq) 

v/|r-r |2 + A 2 
A 



Vlr-roP + A 2 



) 



Vlr-roP + A 2 



(10) 



Here A is the radius of the half-skyrmion at i"o = (xo, yo)- 
Inside the core |r — r | 2 < A, the spin is polarized; out- 
side it |r — ro| 2 > A, one gets a vortex-like spin config- 
uration on XY plane. To stabilize topological vortices 
(half-skyrmions), a small easy-plane anisotropic term 
should be added to the original model phenomenally, 



w 



k J2i n i («>o, f <i) 



From the exact solutions, there exist two types of 
merons : one has a up-spin polarized core Q = — i, 



( n hs)i 



X{x - Xq) 



Kv - vo) 



V|r-ro| 2 + A 2 ' v /| r _ ro |2 + A 2 

A )■ 
V|r-r | 2 + A 2; ' 



. (11) 



+• ~ -r jr It ^ x X 
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FIG. 3: The scheme of (nh s ) 2 , the meron with topological 
charge 1/2. Down-spins locate at the center. 

Also, there are two types of anti-merons : one has a 



up-spin polarized core Q 



( n hs) 3 = ( 



A(y - yo) 



X(x - Xq) 

Vlr-roP + A 2 ' V|r-r | 2 +A 2 

); 



,(13) 



V|r-r | 2 + A 2 ' 



the other has a down-spin polarized core Q 



( n hs) 4 = ( 



X(x - Xq) 



A(y - yo) 



v/|r-r | 2 + A 2 ^/\r-r \ 2 + X 2 

X - ). 

y/\v-r Q \ 2 + X 2 ' 



(14) 
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In AF ordered state, the mass of half-skyrmion vrih, 
associated with the ordered staggered moment : 



m h . 



2np s + E cmc (g < g c 



(15) 



(1 — ■£-) is the renormalized spin stiffness 



where p s = p 

and -Ecorc is the core energy. For the isotropic case n = 0, 
the energy of the half-skyrmion does not depend on its 
radius A and the core energy is zero, -E C orc = 0. For the 
anisotropy case k > 0, the scale invariance is broken at 
H' = k J2 t n i ~ ^kinetic where Akinetic is the kinetic en- 
ergy, Akinetic = |(Vn) 2 . A new scale I 2 = §^ appears 
(a is lattice constant). Now, the core energy is estimated 



as E c 



const 



A\2 



(f) 



which indicates the instability 



of the static soliton against collapse, A — >• 0. The half- 
skyrmion without a core A — > turns into a spin- vortex 
as 



nhs 



x 



A- 



± 



y-yo 



A- 



(16) 

See Fig. 4, all spins of the spin- vortex are suppressed onto 
the XY plane. On the other hand, for the easy-axis 
anisotropic case, k < 0, the energy of spin- vortex will 
diverge, m/ ls ~ Nk — > oo. 



where Vi = ^lTo = ( ipfiA, V'fiB, V^ia, V'J.ib ) and 
V>2 = V4to = ( 4>^2B, ifaiA, $12B, i>\.2A )■ 7m is defined 

1 
1 

t x , t v , t z are Pauli matrices. m e denoted the mass gap 
of the electrons. We have set the Fermi velocity to be 
unit vf = 1. The solution of zero modes has given by 



as 70 = cr C 



<(r) = 



exp(- 
exp(- 

Vo 



k-r | \ 
m e ' 
|r-ro| -i 



and '02 ( r ) 



l£-£o] 



/ 



-exp( 
4 ex P (-^) 

Vo 



■ 56,61 



in Refii 

On the other hand, the double zero modes of a spin- 
vortex on a 21 x 21 lattice was shown in Fig. 5. The ex- 
act charge is localized around the defect center within 
a length-scale ~ A" 1 . For the two spin- vortices, the 
fermion zero modes are slightly split due to tunneling 
between them. From numerical results, we find that 
there exist two zero modes on each spin-vortex. The 
fermion zero modes around a n spin-vortex match the 
results in£&£i. 
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FIG. 4: The scheme of the spin-vortex. All spins are sup- 
pressed onto the XY plane. 




B. Fermionic zero modes 

In this part we calculate the fermionic zero modes 
around the spin-vortex. By replacing the operator Ci and 
site number i by Grassmann number ip(x) and continuum 
coordinates x, y, we calculate the fermion bound state on 
the spin- vortex by the continuum formula of the following 
effective Hamiltonian 

H c g = - J2(Ujclcj+h.c.)+J2 Acjrihs . a c, (17) 

{ij) i 

In the continuum limit, the two-flavor Dirac-like effec- 
tive Lagrangian2i££~— describes the low energy fermionic 
excitations at two nodes ki = k2 = (§, — § ), 

£ C ff = i^iTfi^fi ^1+^27/1^ V>2 +m e (^iiihB ■ (Tipi—fan-hs ■ cr 

(18) 



FIG. 5: The fermionic zero mode of spin- vortex on a 21 x 21 
lattice. 



C. Induced quantum numbers and statistics 

Next, we calculate the induced quantum number on 
the spin- vortices and show the relationship between the 
induced quantum numbers and the topological charge of 
themS 2 -. 

For the solutions of zero modes, there are four zero- 
energy soliton states around a spin- vortex which are de- 
noted by 

) | ai)® | &i>, | a )«) | &o), (19) 

| a )® | 6i), | ai)<8 | & ). 



5 



Here | do) and | foo) are empty states of the zero modes 
?Pi(r) and V'siMi | a\) and | 61) are occupied states of 
them. Without doping, the soliton states of a spin- vortex 
I sol) are denoted by | ao)<8> | &i) and | a\)® \ bo). 

In Refi^ 6 -, the induced quantum numbers on the soli- 
tons states including total induced fermion number Np = 

c\a z Ci and the induced staggered spin number ^ = 

a, i 

\ S c|u z Ci — g c\(J z Ci has been calculated. The to- 

tal induced fermion number on the solitons is zero from 
the cancelation effect between two nodals Np | sol) = 0. 
However, there exists an induced staggered spin moment 
on the soliton states^ 6 -, 5? ^ | sol) = ±i | sol). 

In particular, we will show the relationship between 
the induced staggered spin moment Sf^ ^ and the topo- 
logical charge of the spin-vortex Q. Remember we have 
study the zero modes and the induced quantum number 
on a spin-vortex. Now the spin polarization inside the 
core of the spin- vortex will be determined by the induced 
spin moment : for the case of S 1 ? » = 5 , one has a up- 

spin polarized core, Q= — |; for the case of S? s = —5, 

one has a down-spin polarized core, Q=\- As a result, the 
spin-vortex with induced staggered spin moment becomes 
a half-skyrmion with a narrow core and the topological 
charge of such half-skyrmion dependents on the induced 
spin moment (See Fig. 6). 



In the long range AF ordered state (g < g c ), the spin- 
vortex will be confined and the total energy of it diverges. 
On the other hand, in the region of g > g c , the spin- 
vorterx has finite energy that is E COTC = n. Thus in 
the quantum non-magnetic insulator state {g > g c ), the 
mass of half-skyrmion vanishes in the isotropic limilj^, 
TUhs = and survives when there is a small easy-plane 
anisotropic term mhs = E corc = k. For the easy-axis 
anisotropic case, k < 0, the energy of spin- vortex always 
diverge, m/ ls —t 00. 



IV. QUANTUM SPIN LIQUID STATE 

Whether the non-magnetic insulator is a VBC state? 
Because the VBC state is characterized by the condensa- 
tion of the half-skyrmions, if half-skyrmions are bosons, 
they will condense. Then one gets a VBC state (or quan- 
tum dimer state) 63 ' 64 . However, in the non-magnetic in- 
sulator half-skyrmion obeys fcrmionic statistics by trap- 
ping a bosonic spinon z onto its core. The massless 
fcrmionic vortices lead to a new story of the quantum spin 
liquid states. Consequently, a new type of quantum spin 
liquid state - nodal spin liquid (NSL) is explored. In ad- 
dition, after adding a small easy-plane anisotropic term, 
the fcrmionic vortices get energy gap, then the ground 
state becomes a topologically ordered spin liquid. 



FIG. 6: The spin-vortex with trapped spin moment turns into 
a half skyrmion. The trapped spin moment is denoted by the 
red arrow. 

On the other hand, since each bosonic spinon z carries 
h staggered spin moment, an induced staggered spin mo- 
ment corresponds to a trapped bosonic spinon z. So there 
is a mutual scmion statistics between "bosonic spinon" 
z and the half skyrmion with narrow core (meron or an- 
timeron ). Due to the mutual semion statistics, by bind- 
ing the trapped bosonic spinon z, a mobile half-skyrmion 
becomes a composite fcrmionic particle. 

In order to giving a clear comparation we show these 
quantum numbers of a spin vortex (half-skyrmion with 
narrow core) in Table. 1. 





1 ai)® | b ) 


I a )® | 61) 


Total fermion number 








staggered spin number 


1 

r 1 


1 

\ 


Topological charge 


2 


2 



A. Effective model of half-skyrmions 

Because the spin-vortices may have zero energy, to 
learn the low energy physics of non-magnetic insulator, 
we need to know quantum dynamics of them. In the fol- 
lowing parts we will consider the spin-vortex as a quan- 
tum object and obtain its effective model. 

As shown in Fig. 6, the center of a spin- vortex is the 
plaquette rather than the original sites. So we may define 
dual lattice by /. The spin-vortex will hop from one dual 
lattice to another. Thus the spin-vortices show similar 
behavior of vortices in XY model : it can move on dual 
lattice with the same lattice constant and feel an effective 
7r-fiux phas o 51 ' 52 (See detail in Appendix C). 

We may use the operator to describe such neutral 
fcrmionic particle with half spin at dual lattice /. Here a 
is the spin index. The relation between the zero energy 
states and the fermionic states is given as 

I oi>® I b ) = fl I 0)/ (20) 

and 

I ao)® I &i> = /] t I 0)/ (21) 

(The state |0) / is defined through / Jf |0) f = / n |0) / = 0). 
We call such neutral object at dual lattice I (fermion with 
±i spin degree freedom) a "fcrmionic spinon". Then the 
leading order of the hopping term of the half-skyrmions 



6 



is 

H hs = -J2 (t!,jf}fj + h.a) +J2 m ^(- 1 ) I fifi 

(i,J) i 

(22) 

where // = (fi-f, fii) T are defined as fermionic spinon's 
annihilation operators. / and J denote two nearest- 
neighbor dual sites. Here tj j is the effective hop- 
ping of the fermionic spinons. For the 7r-flux phase for 
the half-skyrmions, one also needs to divide the dual- 
lattice into two sublattices, A and B. After transform- 
ing the hopping term into momentum space, we ob- 
tain a d-wave like dispersion of the fermionic spinons 

6k = ± ^At 2 (cos 2 p x + cos 2 p y ) . So there also exist two 

nodal fermi-points at pi = (§,§), P2 = (§,—§) and 
the spectrum of fermions spinons becomes linear in the 
vicinity of the two nodal points. 

In the continuum limit, we get a Dirac-likc effective La- 
grangian that describes the low energy fermionic spinons 

C hs = ^al^d^a + im hs *a*a (23) 

a a 

where §1 = ^\j = ( foiA, IfiB, /4.1A, /4.1s ) and 
*2 = l4 7o = ( /t2B, ff2A, /i2B, f\.2A )■ In is defined 

(I 0\ 

as 70 = (Jq®t z , 71 = a ®r y , 72 = a ®T x , a = I Q 1 I . 

t x , t v , r z arc Pauli matrices. For simplicity, we set the 
Fermi velocity to be unit. Thus in the isotropic limit 
hi — y 0, we may ignore the mass gap of ^ a and consider 
the fermionic spinons as gapless particles. 



B. Mutual Chern-Simons theory 

From above results, there exist two types of fields, the 
bosonic spinon z and the fermionic spinon /]. Now the 
low energy effective model becomes 

C cS = C hs +C s = — {d^nf+i ^ * QCT (7^ + m hs ) * a(T 

aa 

(24) 

However, there exists non-trivial topological relationship 
between z and $ ai7 - the fields f Q1T that carry ±i winding 
number of AF vector n, 

— e^ A n • x 9 A n = -i V* a — 7 M * a . (25) 

07T ' £ — ' Z 



Here the operator ^- means that the topological charge of 
the fermionic spinons is spin-dependence : For the meron 
with i (— |) spin number, the topological charge Q= — ^ 
(Q=i) with up-spin (down-spin) polarized in the core. 

To ensure such constraint in Eq.(|25j). we add a new 
term in the effective Lagrangian, 

^constraint = ^(^-^An-^nX 9 A n + J V *a7-7M*«)- 
07T ' £ — ' Z 

a 

(26) 

Finally the low energy effective theory is obtained as 



C e tt = Chs + C s + C 



constraint 



^2^a ( 7^ + * y7 M ^ + ^hs) *a + - «a M )z| 2 +C M cs- 



r 



In particular, there exists a mutual-Chcrn-Simons term, 
Cues = -e^ x A^a x = l -^ a »d v A x . (28) 

TT TT 

To obtain above low energy effective theory, we have use 
the equation d v a\ — d\a u = \n ■ d v w x d\\v. 

Then the effective model becomes an U(l) x U(l) 
mutual-Chcrn-Simons (MCS) gauge theory. Fermions 
^aa couple to an U{1) gauge field A M ; boson spinon 
z couples to an <7(1) gauge field a M . The results say 
that act as half topological vortices for a M and z 
act as half topological vortices of A^. This effective La- 
grangian proposed in here and earlier paper o 65 ' 66 retains 



the full symmetries of translation, parity, time-reversal, 
and global spin rotation, in contrast to the conventional 
Chern-Simons theories where the second and third sym- 
metries are usually broken. In the following parts of the 
paper, we use U (1) x U(l) MCS theory to learn the quan- 
tum non-magnetic state near Mott transition. 



C. AF order 

If g < g c , the ground state has a long range Neel order, 
(z) = Zo 7^ or (n) 7^ 0. Hence the effective model turns 
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into 



C cS ~ ^ | a M | 3 + V^a^A 



+i 



where e 2 , ~ 37rm, 



1 

4^ 



and p 2 



p^ is the mo- 



mentum. Are bosonic spinons and fermionic vortices real 

(29) 

quasi-particle? To answer the question, we need to cal- 
culate the energy gap of the gauge fields furthermore. 

Firstly, after integrating we obtain a mass term for 
gauge field Ap, 



where g att is renormalized coupling constant as g g{f = 
■ Here the mass term of the gauge field a M is caused 
by z condensation. After integrating a M we get a renor- 
malized kinetic term for gauge field A^, 



1 

4e A 



1 1 S 



(30) 



with 



4- = ~r^2 • Due to the instanton effect, the gauge 

e A Z Q 7T 

field A^ also obtains a mass gap and a linear confinement 
appears for fermions^. The low energy physics is dom- 
inated only by spin waves and the effective Lagrangian 
becomes £ e ff = 2a~ (<^ n ) 2 - 

"off 

At small T, the AF order is believed in the so-called 
renormalized classical (RC) region. In the RC region, 
the effective Lagrangian loses the Lorentz invariance and 
becomes £ e ff — > ^Ps(Vn) 2 . Here p s is the renormalized 
spin stiffness as p s = p s (l — y). Then at low temper- 
atures, fermionic vortices are paired by the logarithmic- 
attractive interaction V(r) = 2tt/5 s In -[^j . With the in- 
crease of temperature, neutral vortex-antivortex pairs 
like those in the XY model are thermally excited, leading 
to a conventional contribution to the screening effect. By 
the conventional Kosterlitz-Thoulcss (KT) theory, there 
exists a "deconfining" temperature, Td c ~ ^f-- Above 
Tde, free excited fermionic vortices exist. 



z ^"7mA 



1 



{dvArf 



(34) 

Due to exchanging the gauge field A^, a short range in- 
teraction is induced between fermions. It is obvious that 
the short range interaction is irrelevant. As a result, the 
massless Dirac particles \l/ a that couple to A^ become 
real low energy degrees of freedom. 

Secondly, after integrating A^ we get the effective 
model of a M as the 



off 



1 



rr2 V 



(35) 



:(CU M ) 2 

This Lagrangian has the schematic form — 2^2 p 2 + -^P- 
That means the momentum of the gauge field a M is not 
zero, p = ^ 0. Then the gauge field a M has a finite 

2e 2 

energy gap (£ c ff ~ ifra 2 ) and shows roton-like behavior. 
As a result, the induced interaction by exchanging the 
gauge field a M is irrelevant and the bosonic spinons are 
real quasi-particles. 

Finally, the low energy effective theory of nodal spin 
liquid becomes 



C 



cir 



kg ^ 



m 2 z z 2 ] (36) 



1 

4^ 2 



V 2 



D. Nodal spin liquid - isotropic case k = 

If g > <7 C , the ground state has no AF long range order. 
Now z has a mass gap m z , or one has the massive spin-1 
quanta, 



There are three types of quasi-particles : two flavor gap- 
less fermionic spinons ^ a <j, gapped bosonic spinons z 
and the roton-like gauge fields a M . Therefore, from the 
effective theory we conclude that NSL state is stable and 
can be considered as a new type of quantum spin liquid. 



C s = ±- 
2g 



|(d M -*a„)z| +m\ 



(31) 



E. Topological spin liquid - anisotropic case k > 



For the isotropic case k = 0, the fermionic vortex has no 
energy gap. Thus the effective model becomes 



For the anisotropic case k ^ 0, the fermionic vortex 
has finite energy gap. Thus the effective model becomes 



£cff = i^2^a(l^d ll + i^---f fl A t A^ a 



(32) 



ia A1 )z| z + m 2 z 2 



2.9 



+ -e^ x A^a x . 



After integrating the massless fermions and the massive 
bosonic spinons, the effective model turns into 

Ces = -^(d^) 2 ~^ T (d v A^) 2 + -e^ x a^A x (33) 
4e a 4e^ 7r 



(37) 



hi* 9 " 



29 



ia„)z| 2 + m?z 2 



+ -e»» x A u d„a x . 



After integrating the bosonic spinons z and fermionic 
spinons ty a , there appear the kinetic terms for gauge field 
a M and A^, 

C eS = --^ T {d v A, 1 f-^{d li a v f+-e^ x A^d v a x (38) 
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where e 2 



?>irm z and e 2 A ~ §7rm/j S . 
In particular, one may find that due to the mutual CS 
term, the mass term of one gauge field can be obtained 
by integrating the other : by integrating a^, the effective 
model of becomes 



and 



1 e 2 



On the other hand, by integrating we get the effective 



model of as 



-off 



4c 



-(<%a M ) 2 



rr2 V 



Thus the gauge fields have mass gap 



m a = m A 



2tt 



Thus we find that the mutual U(l)xU(l) CS theory de- 
scribes a topological ordered spin liquid without gapless 
excitations. Due to exchanging the gauge fields A^ and 
, a short range interaction is induced between fermionic 
spinons and bosonic spinons. It is obvious that the short 
range interaction is irrelevant. As a result, the spinons 
that couple to A^ or a M become real low energy degrees 
of freedom. 

To show the topological properties of the topological 
spin liquid, we calculate ground state degeneracy on a 
torus. For periodic boundary condition, we can expand 
the gauge fields as 

(A x ,A y ) = (±Q X + J2 Ale*** ±Q y + £ Ale™% 



Pe* 



2tt 



P6y = rUyVy 



2tt' 



The masses are given as M x = -p-jf-, My = ■p-jf- and 
mi = t-t^-, m v = \k^. Thus we map the original 
mutual U(l)xU(l) CS theory to a quantum mechanics 
model of two particles in two dimensions. The effec- 
tive Hamiltonian of (Q x , & y ) and (0 X , 9 y ) corresponds to 
two particles on a torus through two-unit flux : (Q x , 6 y ) 
are the coordinates of the first particle, and (Q y ,6 x ) are 
the coordinates of the second particle. The degeneracy 
for (<d x ,8 y ) degrees of freedom and the degeneracy for 
(6 y , 9 X ) degrees of freedom are given as f (e B ,e v ) = 2 and 
D(& y fi x ) = 2. As a result, for the mutual U(l)xU(l) CS 
theory, the ground states have four-fold degeneracy. That 
means the ground state of the anisotropic case k > is 
a Z2 topological spin liquid. 

For the easy- axis anisotropic case, k < 0, the situation 
changes. For this case, the effective nonlinear a model is 
not available for the temperature below the energy scale 
of the anisotropic k. Thus the ground state are always 
long range (Ising) AF order in the insulator phase M/0. 



F. Experimental predictions 

Firstly, we discuss the experimental predictions in the 
NSL. 

In NSL, the spin-correlation decays exponentially 



a y ) = {j-9 x + J2 «k^- k ±-B y + £ <e- k ) 



(39) 
(40) 



y i where n XtV 



where k ={k Xl k y ) = (-j 

tegers. and (a^a^) are ^ ne § au § e fields with 

non-zero momentum and (Q x ,Q y ) and (6 x ,9 y ) are the 
zero modes with zero momentum for the gauge fields Ai 
and ai. Because the existence of the mass gap, the de- 
gree freedoms for gauge fields with non-zero momentum 
(A£, A^) and (a£, aj£) have nothing to do with the low 
energy physics. The low energy physics is determined by 
(e*,e tf ) and (6 x ,9 y ). 

For the temporal gauge Aq = 0, after the mode expan- 
sion, we write down the following effective Hamiltonian to 
describe the low energy physics of the mutual U(l)xU(l) 
CS theory 



2M X 



2m t 



2m x 



P 2 
2M V 



where the conjugate momentum for (0?, & y ) and (6 X , 9 y ) 
are defined as 



(S + (x,y)S-(0)) 



JQ R Z I + 



(n + (x,y)n-(0)) 



(41) 



with Q = (it, it), r — \J x 2 + y 2 and = n x ± in y . Here 



£ is spin correlated length, £ 



8tt 
Ac/off 



In contrast, in al- 



gebraic spin liquid or algebraic vortex liquid, the spin- 
correlation shows critical behavior. 

Secondly, we calculate the special heat. Because the 
fermionic spinons have no energy gap, the special heat is 
dominated by them. So at low temperature, the special 
heat is 



C v = 



^k 2 T 2 



(42) 



where £ = 



e*+l 



r(3)C(3) ~ 1.803. 



Thirdly, we calculate the spin susceptibility. The defi- 
nition of the spin susceptibility is 



F = F{B = Q)- l - X B 2 . 



(43) 



where F is free energy and B is the external magnetic 
field. There are two contributions to the total spin sus- 
ceptibility x, one from bosonic spinons, the other is from 
the fermionic spinons as 



P e . =M X Q X + 



2tt' 



Pe 



MyQy 



2tt 



X = Xb + Xi- 



(44) 
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Here Xb is given bysi' 



Ki.fiS 



Xb 



2 (2fx B M) 2 . m s p 



7T/3 



ln(e 



1)] 

(45) 

with j3 = 1/fceT. The contribution from the fermionic 
spinons is almost linear temperature dependence as 

41n2 



Xf 



tt;3 



(46) 



in unit of (gfi^) 2 . 

On the other hand, for TSL, the spin-correlation 
also decays exponentially (S + (x, y)S~(0)) ~ e tC ^' R< e _r ^ 
with £ = . However, due to mass gap, rtihs 7^ 0, at 
low temperature, the special heat and the spin suscep- 
tibility from the fermionic spinons are all proportion to 



-pm hs 



and disappear at zero temperature. 



V. CONCLUSION 

Let us draw a conclusion. In this paper, we study the 
non-magnetic insulator state near Mott transition of 2D 
TT-flux Hubbard model on square lattice and find that 
for the isotropic case such non-magnetic insulator state 
is quantum spin liquid state with nodal fermionic exci- 
tations - NSL; for the anisotropy case it is TSL with 
full gapped excitations. The low energy physics is basi- 
cally determined by its U(l)x U(l) mutual Chern-Simons 
gauge theory. There exist both fermionic spinons and 
bosonic spinons. And it is just the mutual semion statis- 
tics between fermionic spinons and bosonic spinons that 
guarantee the stability of quantum spin liquid states. 

Because NSL state represents a new class of quantum 
state which may be applied to learn the nature of the spin 
liquid state in other systems, for example, the Hubbard 
model on honeycomb lattice. People have proposed that 
quantum non-magnetic state of the Hubbard model on 
honeycomb lattice is really a Zi topological spin liquid 
ordered state, which is robust against arbitrary pertur- 
bation including the anisotropic termor—. However, our 
results of the quantum spin liquid state near the Mott 
transition of the 7r-flux Hubbard model arc different - 
for the isotropic case, the non-magnetic state with SU(2) 
spin rotation symmetry is nodal spin liquid (NSL) with 
gapless fermionic excitations and roton-like excitations; 
for the anisotropic case by adding arbitrary small easy- 
plane anisotropic term, the non-magnetic state becomes 
a Zi topological spin liquid with topological degenerate 
ground states. People may check the different theories by 
QMC approach in the future. 

Finally, we give a comparison on different quantum or- 
ders with 7r-vortex (half-skyrmion or vison) and bosonic 
spinon. In general, for a system with 7r-vortex and 
bosonic spinon, there exist five types of quantum orders 

oo69 . 



is always VBC state with spontaneous translation 
symmetry breaking; 

2. AF order: If one has massive bosonic 7r-vortices 
and gapless bosonic spinons, then the ground state 
is an SDW order with spontaneous spin rotation 
symmetry breaking; 

3. Algebraic vortex liquid: if one has gapless fermionic 
7r- vortices and massive bosonic spinons, then 
ground state may be an algebraic vortex liquid 
(AVL)2i. In algebraic vortex liquid, fermionic ex- 
citations themselves couple to massless U(l) gauge 
fields, as gives an example to algebraic spin liquid. 

4. Topological spin liquid: if one has massive fermionic 
7r-vortex and massive bosonic spinon, the ground 
state is a topological order with topologically de- 
generate ground state; 

5. Nodal spin liquid: if one has gapless fermionic tt- 
vortices (the fermionic spinons in this paper) and 
massive bosonic spinons, the ground state is a 
nodal spin liquid without any spontaneous sym- 
metry breaking. In particular, there exist gapped 
roton-like gauge modes. 

In order to giving a clear comparison, we give Table. 2 





7r-vortex 


bosonic spinon 


VBC state 


massless Boson 


massive Boson 


AF order 


massive Boson 


massless Boson 


AVL state 


massless fermion 


massive Boson 


TSL order 


massive Boson(fermion) 


massive Boson 


NSL state 


massless fermion 


massive Boson 



We give a short remark on the relation between NSL 
state and the AVL in Re fi 16 ' 17 i 71 . In AVL, fermionic ex- 
citations coupling to a massless U(l) gauge field and 
cannot be real quasi-particles. In NSL state, due to 
the protection from the mutual semion statistics between 
fermionic spinons and bosonic spinons, fermionic excita- 
tions are real excitations. There is neither bosonic spinon 
nor roton-like gauge mode in AVL state and the low en- 
ergy effective theory of AVL state is also difference from 
that of NSL. In AVL, the spin-correlation shows critical 
behavior, while in NSL state, although there exist gapless 
fermionic spinons, the spin- correlation decays exponen- 
tially. 

This research is supported by SRFDP, NFSC Grant no. 
10874017 and National Basic Research Program of China 
(973 Program) under the grant No. 2011CB921803. 



VI. APPENDIX A: EFFECTIVE NLcrM 



1. VBC state: If one has gapless bosonic tt- vortices 
and massive bosonic spinons, then the ground state 



In this appendix we use the path-integral formulation 
of electrons with spin rotation symmetry to obtain the 
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effective NLcrM of the spin fluctuations^^iL^. The 
interaction term in Eq.[l] can be handled by using the 
SU(2) invariant Hubbard- Stratonovich decomposition in 
the arbitrary on-site unit vector fli 



-{Sli-ctaci} 2 , 



(47) 



Here a = (a x , a y ,a z ) are the Pauli matrices. By replacing 
the electronic operators cj and Cj by Grassmann variables 
c* and Cj , the effective Lagrangian of the 2D generalized 
Hubbard model at half filling is obtained: 

C c e = ^2 c*d T Ci - ^2 (ti,jC*Cj + h.c.) - A ^2 Cj *V<7Cj. 

» (v) i 

(48) 

To describe the spin fluctuations, we use the Haldane's 
mappin g 36 ' 37 ' 39 : 



Sli = (-l)'iWl-L? + L 



(49) 



where = (nf,nf,n*) is the Neel vector that corre- 
sponds to the long-wavelength part of S~2.j with a restric- 
tion nf = 1. L, is the transverse canting field that cor- 
responds to the short-wavelength parts of fli with a re- 
striction Li ■ iii = 0. We then rotate fij to z-axis for the 
spin indices of the electrons at i-site ^ 33 ' 38-41 ' 56 



Ujiii ■ aU l = cr z 
U}Li ■ aUi = h ■ 



(50) 



where Ui €SU(2)/U(1). 

One then can derive the following effective Lagrangian 
after such spin transformation: 



£ ff = ^ 4>*d T ipi + ^2 a o W fa 

i i 

- ^2(Uj^e* a *^ j +h.c) 

<ij> 



%l>i (51) 



where the auxiliary gauge fields = di^\a x + 2a y 
and ao (i) = ao,i (i) u x + ao,2 (i) <J V are defined as 



e ia y =u }u ao(i) = u}d T Ui. 



(52) 



In terms of the mean field result M = (—1)' (i/j*a z il)i) 
as well as the approximations, 



1 - L? ~ 1 



1? 



we obtain the effective Hamiltonian as: 

£ c ff ~ ^2 ^i d Tlpi + ^2 t 00 (*) ~~ A(T ' 

i i 

- /~][U,ji>i(l + ia t j)ipj + h.c] 



aE(-i)V*^ + am^| 



(53) 



By integrating out the fermion fields ip* and ipi, the ef- 
fective action with the quadric terms of [oq (i) — Ac • lj] 



and dij becomes 



2A' 



1?]- 



i r' 9 

5 eff = - / drVh^fao (i)-A(T • l l ) 2 +4 i o s 
•'o i 

(54) 

To give p s and c, for calculation, we choose f/j to be 

^ = ( Z *J f ) , (55) 

where rij = Zjcrz.;, z ; = Zjj,) T , z,Zj = 1. And the spin 
fluctuations around = Zj are 



rij = 2i+Re (00 x+Im (fa) y 



&/2 



Ol 



(56) 
(57) 



Then the quantities UjUj and Ujd T Ui can be expanded 
in the power of fa — cj>j and d T fa, 



UjUj = e-' 1 ^—^ 
Uld T Ui=( t ° 



(58) 
(59) 



\d T fa 

According to Eq. ([52j) . the gauge field and ao (i) are 
given as 

1 

2 



•% - faj) &y 



°0 («) = -ypT&Vy 



(60) 
(61) 



Supposing dij and ao (i) to be a constant in space and 
denoting difa— a and d T fa — iB y , we have 

<kj = -o a " (i-j)^ ( 62 ) 



1 



--B„cr„ 



(63) 



The energy of Hamiltonian of Eq. ([54")) becomes 

S(B y ,a) = ic^ + ^ s a 2 . (64) 

Then one could get ^ and p s from the following equations 
by calculating the partial derivative of the energy 



c 



Ps = 



1 d 2 E (B y) 

N dB 2 y 

1 d 2 E ( a) 
A da 2 



a=0- 



(65) 
(66) 
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k 



Here Eq (B y ) and Eq (a) are the energy of the lower Hub- 
bard band 



(67) 
(68) 



where E^_ k , E_ k and E p + k , E p _ k are the energies of the 
following Hamiltonian US and H p 



Next, to learn the properties of the low energy physics, 
we study the continuum theory of the effective action 
in Eq. (|54[) . In the continuum limit, we denote n^, lj, 
idij ~ UjUj - 1 and a (i) = U}d T Ui by n(x,y), l(x, y), 
Wd x U (or WdyU) and Wd T U, respectively. From the 
relations between Wd^U and 9 M n, 



9 9 9 

a T = a T,l + a T,2 



1 



(d r n) 2 , r = 0, 



2 2,2 

= <Vi + a,, 



■(9 M n) 2 , n = x,y, 



nC = -J2 MWi + h.c) - A^(-1)V> Z ^ 

<ij> i 

+ ^^aoW^ (69) 



ao-l=-- (n x <9 T n) • 1, 



(78) 
(79) 
(80) 



the continuum formulation of the action in Eq. (|54[) turns 
into 



>5 n fF — — 



\J drj A[^ r n) 2 + ft (vnf 



2A 2 



(70) 

Using the Fourier transformations for "H^ , we have the 
spectrum of the lower band of US 



4iA? (n x d T n) ■ 1 + (==- - 4A\)1 21 



51) 



4, k 



where the vector ao is defined as ao = (ao.i, ciq,2, 0) . 

Finally we integrate the transverse canting field 1 and 
obtain the effective NLctM as 



A 2 



(71) 



where £k = ±^J At 2 (cos 2 k x + cos 2 k y ). And £ is obtained 



(82) 



with a constraint n 2 = 1. The coupling constant g and 
spin wave velocity c are defined as 2 - ' 33 ' 38-41 ' 56 : 



k 4 



(72) 



X ± Ps 



r 2_Ps 

X 1 



Similarly, using the Fourier transformations for "H p , we 
obtain the spectrum of the lower band of W 



where x 1 ' is the transverse spin susceptibility 



A 2 + \&\ + \<p\ ± 4A2 1*1 " M* -iV) 

(73) 

where <p and $ are defined as 



A 2 



4(16 



A 2 )§ 



)- 1 -2C/]- 



(83) 



(84) 



y> = -t ^ e* -4 cos 
^-i^e^sinf-a-cJ 



(74) 
(75) 

where 6 =a(e x , e y ) and e 2 = e 2 = 1. Using Eq. (|66)) . p s 
is given as 



Then we use the effective NLctM to study the mag- 
netic properties of the insulator state. The Lagrangian 
of NLctM with a constraint (n 2 = 1) by a Lagrange mul- 
tiplier A becomes 



£cff 



1 

2q? 



(d r n) 2 +c 2 (vnr+zA(l-n 2 ) (85) 



where iX = m 2 and m s is the mass gap of the spin fluctu- 
ations. Using the large- N approximation we rescale the 
field n — >ViVn and obtain the saddle-point equation of 
motion as^~— 



^ k 4 (i&r + A 2 )* 



(76) 



(n ) 2 + k B T 



gc 



cj 2 + c 2 q 2 + m 2 



1. 



(86) 



For the 7r-flux Hubbard model, we get the corresponding 
coefficient e 2 as 



2 _ 1.2 



i 2 [cos (2k x ) (A 2 + 8t 2 + 4t 2 cos (2k v )) 
A 2 + 3t 2 + 1 2 cos (4k x )]. 



In Eq. (|86[) . hq is the mean field value of n and u„ 
2-nnk^T \ n = integers. 

From Eq. (|56"|) . we may get the solution of m s as 



(77) 



= 2fcRTsinh 



e s fc B T sinh 



\2k B Tl 



(87) 
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At zero temperature the solutions of no and m s of Eq. (|86[) 
are determined by the coupling constant g. There exists 

a critical point 3c = x ■ ^ or ^ ne casc °^ 3 ^ X ' we S e t 
a non-magnetic insulator with solutions of tiq and m s : 



7i = (1 - -) 1/2 , m s = 
9c 



(88) 



For the case of g > ^ , we get a long range AF order 
with solutions of no and m s : 



. 1 1, 

no = 0, ?n, = 47tc( ) 

3c 3 



(89) 



VII. APPENDIX B : THE INDUCED KINETIC 
TERM OF GAUGE FIELD 



Starting from the CP(1) model, 



1 r 

23 



\{dft - m M )z| + m 2 z^ 



we calculate the induced the kinetic term of the gauge 
field a ft- 

We obtain the expression of the expansion of the renor- 
malization 2-point function at one-loop order : 



(27T) 



(Ay, + 2p M )(fc„ +2p„) 
(k 2 + rn 2 )((fc + p) 2 + rn 2 ) 



-2(5 



1 



(90) 

At low energy limit as p 2 — > 0, we simplify the above 
integral to the following one 

^rlTli^-W 2 ) (9i) 



(27T)' 

2 



2 J 3m 2 
{VpPv - S^p 2 ) 



At low energy limit ( p 2 — > 0). Thus from — ^(c^a^) 2 , 
the induced coupling constants of three dimensional 
gauge field is obtained as 



Z o z 



(92) 



Using similar approach, oen may get the induced ki- 
netic term of gauge field A^. 



VIII. APPENDIX C : PROJECTIVE 
SYMMETRY GROUP OF HALF SKYRMIONS 

In the insulator state (non-magnetic or magnetic or- 
der), there exists a bosonic spinon on each site. Because, 
the half skyrmion is really the vortex, we may consider 
the dynamics of half skyrmions as that of the vortices of 



the insulator state of bosons on lattice with unit filling 
(See Fig. 7). As a result, the square lattice of bosonic 
spinons plays a role of 7r-fiux phase of half skyrmion on 
dual lattice I = (I x , I y ). So we may use the approach of 
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FIG. 7: Spin-vortex on dual lattice 



the projective representation of the space group (PSG) 
for vortex to learn the dynamics of half skyrmion here. 
The PSG denotes the phase transformations associating 
the operations in the space group. 

For the 7r-fiux phase of half skyrmions, we have x and 
y translations, and a 7r/2 rotation of the PSG on the dual 
square lattice as 



T 



R 



tt/2 



/„) 

Iy) 

f{I X , Iy) 



/(4 + 1, Iy)(-l)' v 
f(I x ,-I y )(-lY' +I > 



(93) 



where f(I x ,I y ) is the wave-function of half skyrmions 
and I = (I x ,I y ) are the dual lattices. One may check 
that the operations associated with translations in the 
PSG do not commute as T x T y = —T y T x that means the 
half skyrmion obtains a 7r phase encircling a site of the 
direct lattice (i) with one bosonic spinon. 

We are looking at the half skyrmions hopping around 
dual square lattice in the presence of 7r-flux per plaque- 
tte. Then the effective action should be invariant under 
PSG and the dual U(l) gauge symmetry. For the sim- 
plest gauge and nearest neighbor hopping, the effective 
Hamiltonian of half skyrmions is 



(i, J) 



ii,j/i T /j + h.c. 



(94) 



where the effective nearest-neighbor hopping ti t j could 
be chosen as h,i+x = t, U,i+y = te ±l ^ . 
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